Majorana qubit rotations in microwave cavities 
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Majorana bound states have been proposed as building blocks for qubits on which certain oper- 
ations can be performed in a topologically protected way using braiding. However, the set of these 
protected operations is not sufficient to realize universal quantum computing. We show that the 
electric field in a microwave cavity can induce Rabi oscillations between adjacent Majorana bound 
states. These oscillations can be used to implement an additional single-qubit gate. Supplemented 
with one braiding operation, this gate allows to perform arbitrary single-qubit operations. 
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The search for Majorana bound states (MBS) [T] in 
solid-state systems is currently receiving a lot of attention 
from theorists and experimentalists alike [2H1]- Shortly 
after the prediction that MBS can be realized in semicon- 
ductor nanowires with strong Rashba spin-orbit coupling 
[SHZ], various experiments have indeed reported signa- 
tures of Majorana fermions in such systems, which are 
currently being scrutinized [BlUl|. MBS appear at the 
phase boundaries between topologically trivial and non- 
trivial sections of the nanowire. The latter can be cre- 
ated in the presence of a proximity-induced supercon- 
ducting gap A and a magnetic field B: if /i denotes 
the (position-dependent) chemical potential, the wire is 
in the topologically nontrivial (trivial) phase in regions 
where < B 2 ~ A 2 (p 2 > B 2 - A 2 ) 0. 

MBS are interesting in their own right because they 
constitute the simplest quasiparticles with non-Abelian 
exchange statistics. They have also been suggested as 
resources for topological quantum computing J2j where 
the fact that a single qubit can be encoded in two or four 
spatially separated Majorana fermions offers protection 
from some common sources of decoherence. 

Current proposals for topological quantum computing 
are striving to overcome two obstacles: First, it has been 
shown that the protection of MBS from decoherence in 
realistic systems is less than ideal [T5HT5] . Indeed, quasi- 
particle poisoning is probably a ubiquitous source of de- 
coherence in all proposed setups to date. If one assumes 
that this problem can be overcome, a second, more fun- 
damental difficulty arises: topologically protected opera- 
tions on Majorana fermions, i.e., braiding, are not suffi- 
cient to achieve universal quantum computing 12] ■ They 
have to be complemented by other single-qubit and two- 
qubit gates - as well as read-out operations - which will 
still have to be performed in a topologically unprotected 
way. 

In this paper, we shall focus on the second problem. 
In order to find a way to perform operations on a Majo- 
rana qubit in a topologically unprotected, but minimally 
invasive way, we shall investigate two MBS which are 
coupled via a gapped system of length L. This gener- 
ally leads to an energy splitting between the two MBS 



which is proportional to e~ L ^, where £(/i) is a model- 
dependent decay length which depends on the band gap 
and increases for /.i — > 0. We shall show that by cou- 
pling the MBS to an electric field, it becomes possible to 
tune £, thus allowing for a coherent manipulation of the 
MBS. 

In order to see the benefits of such an endeavor, let us 
consider a popular proposal for a Majorana-based qubit. 
Four Majorana modes 71.2,3,4, which satisfy 7J = jj and 
{ijilk} — 25jfc , can be used to define two Dirac fermion 
operators ip L = (ji + «72)/2 and tp R = (73 + iji)/2. 
Since the topological protection relies on a conserved 
fermion parity, the computational basis should contain 
states with the same parity. Therefore, one can use 

ll)=4|0), |t)=4lo) (1) 

as the two logical states of the qubit. Certain topo- 
logically protected single-qubit gates can be realized by 
braiding. For instance, exchanging the positions of the 
MBS 71 and 72 (73 and 74) corresponds to the transfor- 
mations 

(Z7T \ / ITT \ 

^o- z J, I/34 = cxp \—£<r z J , (2) 

respectively, where a z is a Pauli matrix in the basis 
{It) 1 II)}- Similarly, exchanging 72 and 73 corresponds 
to U23 = exp (— ma x /A). However, it is easy to see that 
these three operations are not sufficient to reach arbitrary 
points on the Bloch sphere spanned by all normalized lin- 
ear combinations of |f) and \\). 

In order to perform arbitrary single-qubit rotations, 
one needs to supplement this set of gates by topologically 
unprotected operations. If the MBS 72 and 73 are subject 
to a coupling Hamiltonian 

is 

V = —7273 (3) 

for a certain time t, it corresponds to the following uni- 
tary transformation of the qubit state 

U e (t)=ew(- t f<rX (4) 
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FIG. 1. (Color online) Four Majorana bound states (MBS, 
stars) 71,2,3,4 at the edges of topological superconductors (TS, 
dark blue) encode one logical qubit. The central MBS 72 and 
73 are tunnel coupled (t c ) to a topologically trivial, gapped 
central region (CR, light yellow) with length L. All energies 
are small compared to the induced gap A. 

Arbitrary points on the Bloch sphere can now be reached 
by combining this operation with braiding. For instance, 
the action of t/ e (ti)i/i 2 i7 e (i 2 ) on the initial state |f) pro- 
duces a final state with components, 

(a x ) = -sin(et 2 ), 

(cr y ) = -sin(eti)cos(et 2 ), 

(a z ) = cos(eii) cos(rf 2 )- (5) 

The vector ((u x ) ,(<J y ) ,(<t 2 )) indeed covers the entire 
Bloch sphere for et 2 € [-tt/2, tt/2] and eh € [0, 2vr]. The 
same is true for all initial states on the Bloch sphere, so 
one can realize universal singlc-qubit rotations using the 
operator U e (t) complemented by the braiding operation 
Ui2- Therefore, our aim is to find a physical mechanism 
to realize the coupling Hamiltonian Q in a controllable 
way. 

Our starting point is the minimal system which con- 
tains the two MBSs 72 and 73, which are coupled indi- 
rectly by tunneling to a common, gapped central region 
(CR). The system is depicted schematically in Fig. [I] It 
can be described by the Hamiltonian 

H Q = J2 <k)d\d k 

k 

| 72 [<i(0) - rft(0)] - l f[d(L) + d\L)] lz . (6) 

We assume that the generic single-particle spectrum of 
the CR is e(fc) = t 2 /(2m)-fi, where [i < 0. Here, k is the 
wave number which labels the eigenmodes dk , and m > 
is an effective mass determined by the band curvature. 
The CR has a length L, and tunneling (with amplitude 
t c ) from and to the MBS occurs at x = and x = L. 
In order to impose the proper boundary conditions, we 
use the conventional "unfolding" transformation [TB] to 
map the CR onto a chiral system with length 2L con- 
taining only right-movers. One can then impose peri- 
odic boundary conditions on this doubled system, which 
leads to the momentum quantization fc = 7rn/L, where 
n > 0. The fermion operators in real space are given by 
d(x) = (l/V2lm k e ik *d k . 



The Hamiltonian Hq is quadratic and can easily be 
solved exactly. One finds that tunneling leads to a 
nonzero overlap between 72 and 73. For small t c , the 
resulting level splitting due to virtual cotunneling pro- 
cesses is 



with the decay length £ = (2m|/i|) -1 / 2 . At energies small 
compared to the effective coupling between the MBS 
can be described by the Hamiltonian ([3]). Therefore, the 
operation U s (t) could in principle be realized by tuning 
the chemical potential /1 time-dependently. However, it is 
unlikely that this can be done on short timescales and ac- 
curately enough to realize qubit rotations reliably. More- 
over, tuning I /1 1 to small values is undesirable. Once real 
tunneling processes become possible (/1 > 0), the Majo- 
rana qubit will quickly decohere [14 . 

A possible way to realize a coupling between MBS in 
a more controllable fashion is to use microwave photons 
to activate the tunnel process. Indeed, the induced gaps 
in recent experiments are in the microwave regime [8]. 
It may at first be surprising that Majorana fermions, 
which are by definition uncharged, should be suscepti- 
ble to electromagnetic radiation. However, today's solid- 
state versions are superpositions of particles and holes, 
and are chargeless only on average. The existence of a 
coupling of MBS to a vector potential A(x) can be de- 
duced explicitly by using the minimal coupling substi- 
tution p — > p — eA{x) in semiconductor nanowire based 
proposals [6j [7] , or equivalently by using the Peierls sub- 
stitution t — > te~ ie J dxA ( x ) i n a Kitaev chain [T7j, which 
constitutes the simplest exactly solvable model hosting 
MBS [T]. 

Therefore, we place the previous system in a microwave 
cavity and describe it using the Hamiltonian, 

H = <k)d\d k -g c (a + a^) 72 [d(0) ~ t*t(0)] 
fc 

- ig c {a + a))[d(L) + d) (L)} l3 + H ph . (8) 

In contrast to Eq. ([6]), a tunnel process between the CR 
and the MBS is now accompanied by the emission or 
absorption of a photon. The cavity degrees of freedom 
are described by the Hamiltonian _ff p h = tta^a + -ffbathj 
where f2 is the resonance frequency and -ffbath accounts 
for driving and damping of the cavity. A similar form of 
activated tunneling into MBS was investigated previously 
in Ref. [18]. 

In order to determine the effective coupling between 
the MBS due to the electron-photon interactions, we cal- 
culate the time-ordered Green's function, 

£> 72 (i) = -i<T 72 (t) 72 (0)), (9) 

where T denotes the time-ordering operator. Since the 
Hamiltonian H is no longer exactly solvable, we use a 
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FIG. 2. Feynman diagram for the lowest-order terms of the 



Dyson equation (j 1 Op . D- h ., G;j, and A denote the Green's 



functions of the Ivlajorana bound states, the central region, 
and the photon, respectively. 



perturbative approach valid at small coupling strength 
g c and perform a resummation of Fock type diagrams, 
see Fig. [2] The result for D l2 (u) is 

D l2 = 4°! + £K°)£ 22 £> 72 + 4^230(^32^, (10) 



where [Zr 1 ^ 
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S33. The unperturbed Ma- 



jorana Green's functions are trivial because for g c = 0, 
72,3 are time- independent, i.e., D^ s (t) = — isgn(t). The 
self-energy is due to tunneling with absorption or emis- 
sion of a photon, Sij(t) = ig^A(t)Gij(t), where A(t) and 
Gij (t) are Green's functions of the photon and the CR, 
respectively, which will be discussed in the following. 
The photon Green's function is defined by 

A(t) = -1 (To(t)ot(0)) Q - i (Tat(*)a(0)) . (11) 

Let us first assume that the cavity field is coupled to 
a thermal environment [19 . The time-dependence of 
a(t) is then governed by the Langevin equation, a = 
(—if) — n/2)a + Y^ain, where k is the cavity decay rate 
and a- ln (t) is the fluctuating thermal mode, satisfying 

(a[ n (w)ain(w')) — 2TTn t h(u>)5(uj — ui'), where nth(^) is the 
Bose distribution. Approximating n t h(w) « nth(^) =: 
rtph, one finds for n p h 3> 1, 



7)=± 



(k/2) 2 + (w - 77ft) 2 



(12) 



The same result formally applies if we assume that the 
cavity is coherently driven at its resonance frequency f2, 
but in that should be interpreted as the line width 

of the input field and n p h is proportional to the number 
of photons in the cavity. For n = 0, A(uj) consists of two 
delta peaks at u = ±0. A nonzero line width turns the 
latter into Lorentzians of width k. 

Finally, G^{€) is the unperturbed Green's function of 
the CR. Defining G d (x,x',t) = -i(Td(x, t)<#{x', 0)) , 
one finds 

G 22 (t) = G 33 (t) = G d (0, 0, t) - G d (0, 0, -t), (13) 
G 23 (t) = -G 32 (-i) = ~iG d (0, L, t) - iG d (L, 0, -t). 



The Fourier transform of the Green's function 
G d (x,x' , t) has the standard form, G d (x,x',uj) = 
2X Efc 0(k)e lk( > x - x "}/[uj - e(fc) + iQ+], where we took into 
account that sgn[e(fc)] = 1 due to the gapped spectrum. 
Equations (10l-((l3| can be used to obtain the MBS 



Green's function, 
■ dto 

D l2 (t) 



27r [I-E22H] + E 23 (u,)£ 23 (- 



(14) 



The residue theorem states that this Fourier transform 
at positive t should be determined by the poles of the 
integrand in the complex lower half plane. For g c = 0, 
there is a single pole at u = 0. For small g c , the pole is 
split and shifted by an amount proportional to g\. Since 
TiijioS) are regular functions near u — 0, we can replace 
Ejj(o>) by £,j(0). These functions are given by E 22 (0) = 
and 



E 23 (0) 



-in ph glmL ^ 
r,=± 



1 



^sin(y / ^ ? ) 



(15) 



where we defined <j>± = (±Q + /i + in/2)/eL and = 
(2mL 2 ) -1 is proportional to the energy of the lowest ex- 
cited state in the CR. The Fourier transform ( 14 ) now 
leads to 



D l2 (t>0) 



_ ie -2\r R \t e -2is g n(r R )n R t 



(16) 



where = Re[E 23 (0)] and Qr = Im[E 23 (0)] denote real 
and imaginary parts of the self-energy, respectively. The 
Green's function thus displays damped Rabi oscillations 
of the quantum state of the MBS. The frequency \Qr\ 
and damping rate \Tr\ are determined by the relation 
between the energy scales e^,, f2, k, and /i. 

First, let us focus on the regime k <C ^ < \fi\. In this 
case, the photon energy is insufficient to overcome the 
band gap. For weak damping, we find 



ir H | = 



\n R \ = 



Up^glnL 

"■phff 2 



exp 



exp 



-Vi-o/|m| 



(17) 



These functions are shown in Fig. [3] The damping rate 
r r is determined by the photon line width k. For fixed k, 
the ratio between damping and Rabi frequency is given 
by IFr/Or] = n/(4y/\Q + /i|cl). The Rabi frequency is, 
as expected, exponentially suppressed in the length of the 
CR. However, as the photon frequency ft approaches the 
critical value the prefactor ^/l — < 1 leads to 

a significant increase of £Ir. As long as Mr <C \/i\, these 
Rabi oscillations lead to an effective low-energy Hamilto- 
nian of the form ^ with e = flR. Therefore, single-qubit 
rotations can be achieved by driving the system for a fi- 
nite time with a microwave frequency close to the gap. 

Next, we discuss the case > In this regime, the 
photon field can excite real electrons and holes from the 
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neling can have a strong impact on the coupling between 
two adjacent Majorana bound states (MBS) which are 
separated by a gapped region. Absorption and emission 
of photons cause Rabi oscillations of the MBS with a fre- 
quency that depends sensitively on the difference between 
the photon frequency and the gap width. Damping of the 
Rabi oscillations is only caused by a nonzero photon line 
width. For subgap photon frequencies ft < the Rabi 
frequency is a monotonic function of ft and increases ex- 
ponentially for ft — ¥ Such Rabi oscillations can be 
useful when applied to Majorana-based qubits because, 
complemented by braiding, they allow for the implemen- 
tation of arbitrary single-qubit rotations. 

This work was financially supported by the Swiss NSF 
and the NCCR Quantum Science and Technology. 



FIG. 3. (Color online) Rabi frequency \Qr\ and damping rate 
\Tr\ in units of c = n^gl/L for /i = — lOOei, and k = 5el. 
A large photon line width re has been chosen to highlight the 
essential features. The crosses are the numerical solution of 
Eq. \1A\ . Solid green and red lines correspond to the solutions 
for the limits of a low- frequenc y p hoton, see Eq. (17 1, and 
high-frequency photon, see Eq. (18 1, respectively. 



MBS into the CR. One still finds £22(0) = 0, but the 
self-energy £23(0) now depends sensitively on the level 
structure of the CR. Since the topological protection of 
the MBS relies on a large length of the CR, let us focus 
on the limit <C -C K <C ft. In this case, the sine 



function in the denominator of Eq. (151 causes Lorentzian 



resonances whenever the photon frequency matches an 
eigenenergy of the CR. For ft n 2 7r 2 e^ — fi, 



|r*| 
in*! 



"ph.9 c 

L (ft- 



phSc 



(«/2) 

- /! - n 2 n 2 e L ) 2 + (k/2Y 
(ft + fi- n 2 TT 2 e L ) 



L (n + fi-n 2 Tr 2 e L ) 2 + (K/2) 2 ' 



(18) 



The resonances are shown in Fig.[3j Finally, let us remark 
that similar Lorentzian resonances in F^j and ft^ as a 
function of ft also arise in the limit of a very short CR, 
where one can replace the CR with a single fermionic 
level at energy — \i > 0. 

Conclusion. We have shown that photon assisted tun- 
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